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Using tlie non-equilibrium Keldysh Green's function formalism, we show that the non-equilibrium 
charge transport in nanoscopic quantum networks takes place via current eigenmodes that possess 
characteristic spatial patterns. We identify the microscopic relation between the current patterns 
and the network's electronic structure and topology and demonstrate that these patterns can be 
selected via gating or constrictions, providing new venues for manipulating charge transport at the 
nanoscale. Finally, decreasing the dephasing time leads to a smooth evolution of the current patterns 
from those of a ballistic quantum network to those of a classical resistor network. 

PACS numbers: 73.63.-b, 73.22.-f 
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Understanding charge transport in nanoscale systems 
[l|-|j| has attracted significant interest over the last few 
years in the context of molecular electronics [5| and 
through recent advances in the fabrication of artificial 
quantum structures |6l4llj . The coherent nature of such 
nanoscopic systems is of particular importance: it leads 
to the formation of eigenmodes in the density of states, 
which are the basis for the engineering of novel quantum 
phenomena, such as quantum imaging [lOj . The question 
naturally arises of how the interplay between the spatial 
structure of such eigenmodes, dephasing, and the leads' 
location and width affects the charge transport through 
nanostructures. 

In this article, we address this question using the non- 
equilibrium Keldysh Green's function formalism [l2|, [l3J. 
We show that nanoscopic quantum networks possess cur- 
rent eigenmodes exhibiting distinct spatial current pat- 
terns that can be selected via gating or through con- 
strictions. We identify the microscopic relation between 
these current patterns, the electronic structure (and ge- 
ometry) of the quantum network, and the leads' location 
and width for different network topologies. Moreover, 
we show that the existence of circulating current loops in 
these eigenmodes can be described via their enstrophy, 
and that modes with a large enstrophy in general carry a 
small total current. Finally, we show that with decreas- 
ing dephasing time, the current patterns evolve smoothly 
from those of a ballistic quantum network to those of clas- 
sical resistor network. These results provide important 
insight into the question of how charge transport can be 
controlled and manipulated at the nanoscale. 

We consider a two-dimensional quantum network con- 
sisting of A^ = Nx 



He 



X Ny sites with Hamiltonian 



r,r', CT 1,1', cr 

- ^ Wofltflr + g ^(4 -t- ar)4 ,jCr,CT 



Here, c]. ^ {^\ cr) creates an electron with spin a at site 
r in the network (1 in the leads), i, t;, and t/j are the 
hopping matrix elements between neighboring sites in 
the network, in the leads, and between the network and 
the leads, respectively. Below, we set for concreteness 
ti/t = 10, and th/t = 0.1. The last term describes the 
interaction of electrons with a local phonon mode with 
energy wq. 

The current between adjacent sites r,r' in the network 
is induced by different chemical potentials, itL,R in the 
left and right leads, and given by [13| 
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Here, G^ is the full Keldysh Green's function matrix 
of the network and leads, which accounts for the elec- 
tronic hopping, the leads' electronic structure, and the 
electron-phonon interaction. For a non-interacting net- 
work (.9 == 0), one has G^ = (l - g^'i)'^ g^ (l - ig")"^ 
with. g^{uj) == 2i[l~2hFiuj)]lm[g'^{u})]. Here, g'''"''^ 
are the retarded, advanced and Keldysh Green's function 
matrices, respectively, containing the Green's functions 
of the decoupled {t, th — 0) network sites and leads. This 
formulation of the problem is sufficiently general such 
that the network sites can represent atoms, molecules or 
quantum dots, the only difference being the form of g. In 
what follows, we will assume for simplicity that each site 
contains only a single, relevant electronic level, such that 
g^ = l/(a; -I- iS) with 5 = 0+, and will compute the lead 
Greens functions using the renormalization procedure of 
Ref. [l_4] . Finally, np is a diagonal matrix containing the 
Fermi-distribution functions (with fcsT — lO^^i unless 
otherwise stated) and t is the symmetric hopping ma- 
trix. 

We begin by studying charge transport through the 
(11 X 21) non-interacting quantum network {g = 0) 
shown in Fig. [TJ We consider the ballistic (quantum) 
limit, where I^r' does not depend on the chemical po- 
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FIG. 1. (a) Conductance, G and the enstrophy, 77. Spatial current patterns for A^ = O.OOlf, 5 = 10~"^t, and (b),(g) fj,c ~ 0, (c) 
fic = O.llSi, (d) fic = 1.819i, (e) /ic = 1.932i, (f) fic ~ 3.912f. (h) Real and (i) imaginary part of G^^l at oj = 0, normalized by 
their maximum value. Spatial current patterns for /ic = in a (j) (11 x 13), (k) (11 x 19), and (1) (11 x 25) network. Currents 
are normalized by Imax ~ max\Iry\. 



tential in the network, and current conservation is au- 
tomatically satisfied [13|. When the network is discon- 
nected from the leads, it possesses discrete states with 
energies -E'(n) ~ — 2t [cos(/i;„j^ao) -I- cos(fc„^ao)] where 
km = niTT/{Ni + l)ao, Ui^ !,■■■ ,Ni {i = x,y), and 
ao is the lattice spacing [J]. To identify how these states 
contribute to the charge transport when two leads are 
connected to single network sites labeled L and R [see 
Fig. mb)] , we plot in Figdja) the conductance of the net- 
work, G{Vc) ~ I{Vc)/ ^V in the limit of vanishing bias 
AV = [jiL — /^fl)/e — )■ 0, as a function of bias midpoint 
Vc — [i'd E where [ic = (ml + /ifl)/2. In this limit, one 

obtains Giy^) = ^'^^tiNl\G\^;R,{lJ^c)? , where iVo is the 
leads' local density of states, and G\ p^ is the non-local 
Green's function between L and R [13]. G{yc) exhibits 
a resonance (i.e., a peak) whenever [i^ coincides with the 
energy of a network's eigenstate, -E(n), whose wave func- 
tion does not vanish at L and R. We refer to these 
resonances of the conductance as the current eigenmodes 
of the network. They are the non-equilibrium analog of 
the eigenmodes in the (equilibrium) density of states [10| , 
with both types of eigenmodes possessing distinct spatial 
patterns [J_5]. Note that for the lead positions shown in 
Fig.[TJb), only 121 states (those with odd 
network's 231 states can carry a current. 



Uy ) out of the 



Each of the current eigenmodes exhibits a distinct spa- 
tial current pattern as shown in Figs.[IJb)-(f) [the current 



patterns shown in Figs.[IJb) and (c) correspond to modes 
(1) and (2) in Fig. [Ija), respectively]. Hence, by gating 
the network (i.e., varying fj.,, ) one can not only alter the 
total current flowing through the network, but also the 
spatial path that it takes [3|, |j, |9| . The current patterns 
in Figs, mb) - (e) exhibit loops of circulating currents 
[SJ, y, which are a characteristic feature of the quan- 
tum (ballistic) limit, and are absent in classical networks 
(see below). Note that the current loops in Fig.IlJb) are 
detached from the main current path, and therefore do 
not contribute to the total current through the network. 
These loops might be experimentally detectable through 
the magnetic dipole fields they generate inside the net- 
work, in particular, when, as is the case in Fig.[TJb), they 
create a dipole field that is opposite to the one generated 
by the total current. Associated with the current loops 
are links in the network in which the current flows op- 
posite to the applied bias [see center row of Fjg. [2b)], 
a phenomenon referred to as current backflow jl6| . Note 
that the current pattern for a finite applied bias AV is 
a superposition of all current patterns associated with 
current eigenmodes lying between /i^ and /i^, as follows 
directly from Eq.(l2]). 



In order to understand the spatial structure of these 
current patterns, we note that for T = and AV -^ 0, 
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FIG. 2. Current patterns for S = 10" i, fic = 0, A/i = 
2 X 10~''t (a) for narrow leads, and (b),(c) wide leads with 
constrictions. 



the current between two neighboring sites, r, r' is 



I te 



/rr' - ANo4 — Im [C^j^Gl,,] 



AV 



(3) 



The contribution from hopping to the right lead is in- 
cluded in the above expression since we made use of the 
identity Im [C^j^Gl^,] = -lm[Gl^G^^,]. This form of 
Jrr' demonstrates the non-local character of charge trans- 
port in a quantum network: the current between sites r 
and r' arises from electrons that first propagate from r to 
the left lead and then to site r' [as indicated by the yellow 
arrows in Fig. [Ifg)], rather than from electrons directly 
hopping between the two sites. Eq.® also reveals the 
microscopic origin of the current's spatial form in more 
detail. In particular, a comparison of the current pat- 
tern shown in Fig. [ijg) with the corresponding real and 
imaginary parts of the non-local Green's function, G^l, 
plotted in Figs. [Ijh) and (i), respectively, demonstrates 
that the current path is primarily determined by Re G^l • 
Therefore, the current pattern of a given eigenmode does 
not reflect the spatial structure of the wave-function of 
the associated eigenstate, which is reflected in the spatial 
form of Im G^l , in contrast to earlier findings ^, ,£] . As a 
result, the spatial forms of the network's local density of 
states, A^(r,a;) = — ImGj; ^ ^/n and of the current eigen- 
mode at uj — fic are in general not related [note that in 



the hmit th —> 0, one has [ImG^ ^j = jlmG^ lI for w = 0, 
and therefore in this limit iV(r, w = 0) possesses the same 
spatial form as Im Gli^{Ld = shown in Fig. [TJi)]. This 
conclusion holds in general for all networks and eigen- 
modes we have considered so far. It is interesting to 
note that the eigenstate contributing to the current flow 
for /ic = [Fig. [TJb) and (g)] possesses the wave-vector 
{kn^,kny) — (7r/2,7r/2), and one might therefore expect 
that the current propagates (on average) along the diag- 
onal direction. However, such a current pattern is only 
realized if the two leads can be connected by a "diag- 
onal" path, which is possible for the case in Fig. [TJg) 
[where the path resembles the ballistic propagation of an 
electron that bounces off the network boundaries], but 
not for the case in Fig. [ijb). Finally, we note that the 
spatial current patterns exhibit a very rich behavior as a 
function of the network's aspect ratio Ny/N^, as shown 
in Figs. [TJb), (j)-(l) for fixed iV^; = 11, and varying Ny 
at ^c — 0. 

To determine the effect of current loops on the total 
current carried by a mode, we consider the normalized 



enstrophy 77(14) 



Er|VxIr(K)|V4fHere, the 



sum runs over all plaquettes of the network, and V x Ir 
is the curl of the current around a plaquette centered at r. 
77 is a measure for the vorticity of the current pattern and 
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FIG. 3. Current patterns on a cylinder (the cyhnder is cut 
along the axis and flattened out to show the currents) with 
He = 0, A/i = 10"^f, and S = lO'^t. 



thus the extent of the circulating current loops. A com- 
parison of T] and G [see Fig. [Ija), for clarity, only a small 
range of Vc is shown] reveals that a current mode with a 
large enstrophy in general carries a small total current. 
In cases where there are more than one state located be- 
tween fj.L and fiR, one needs to consider the enstrophy 
and total current carried by each of these states. Note 
that for the large enstrophy mode shown in Fig.[ljc), the 
currents inside the network are significantly larger than 
the total current. 

To understand the effect of the leads' width on the 
charge transport, we consider the (29 x 9) network shown 
in Fig. [2] at /Ltc = 0. For narrow leads [Fig. [Ha)], the 
current follows a ballistic diagonal path connecting the 
two leads. In contrast, for wide leads the current flows 
almost uniformly along the rows of the network, with 
a larger current flowing along the edge rows [Fig. Hfb)]. 
Such a uniform current pattern can be considered as a 
superposition of current patterns associated with narrow 
leads. As a result, one can select narrow lead current 
patterns [see Fig. [SJa)] in wide lead systems by using 
constrictions as shown in Figs, [^c) and (d), in analogy 
to mesoscopic systems [yj. Moreover, a large degeneracy 
of a current carrying state exerts a subtle, local effect 
on the current pattern, as follows from a comparison of 
Fig. IHg), where the E = eigenstate is non-degenerate, 
and Fig. EJa), where the E = eigenstate is nine- fold 
degenerate. While in both cases, the current flows along a 
diagonal direction, it flows along both sides of a plaquette 
of width ao for the (29 x 9) network [see yellow plaquette 
in Fig.[2l[a)]), while for the (11 x 21) network, it flows only 
along a single side of a plaquette of width 2ao [see yellow 
plaquette in Fig. [Hg)]). This change in the local current 
patterns arises from the different spatial dependence of 
Im [GLr] in both cases. 

To demonstrate that the above qualitative features of 
current flow are robust against changes in the topology 
of a network, we consider in Fig.|3]a cylinder network (as 
are realized by carbon nanotubes [7|) with narrow leads 
at fie = 0. For a cylinder with circumference Ny = 20 
and length N^ = 59 [Fig. ^a)] (where the E = state is 
18-fold degenerate), the two leads are again connected by 
a diagonal current path that winds around the cylinder; 
its local and global spatial patterns are similar to those 
shown in Fig. HJa). When one of the leads is rotated 
(along the circumference of the cylinder) by an angle of 
IT [see Fig. EJb)], it is no longer possible to connect the 
two leads by a diagonal path, and one obtains a super- 
position of two current patterns [that shown in Fig. |3l[a)] 
which differ by a tt rotation around the cylinder's axis. In 
contrast, when the circumference of the cylinder is short- 
ened, as shown in Fig. ^c) for a (59 x 18) network, the 
leads cannot be connected by a diagonal path, and the 
resulting current pattern exhibits a significant amount of 
transverse current flow, i.e., current flow perpendicular to 
the direction of the applied bias, as well as current back- 



flow. Only when one lead is rotated by tt/3 along the 
circumference with respect to the other lead, as shown 
in Fig. IS^d), does the current flow predominantly along 
a diagonal path connecting the two leads, though some 
transverse current flow remains. 



In order to study the effects of dephasing [3, ll7| on the 
spatial current patterns in nanoscopic networks, we con- 
sider the electron-phonon interaction [see Eq.([T])] [18| and 
employ the high-temperature approximation ujq <IC ksT 
(with ujQ — > 0) introduced in Ref. [111]. In this case, 
the self-energy is related to the full Green's function via 
S^''' = iG^''', where 7 = 2g'^T/uJo, such that the de- 
phasing process is controlled by a single parameter, 7 
(a derivation of the Green's functions in this approxi- 
mation is presented in Appendix A). The self-consistent 
solution of the resulting Dyson equations for G^'*" guar- 
antees current conservation in the network. In Fig.|4]we 
present the evolution of the current pattern for a (7 x 13) 
network with increasing 7. As 7 increases, the current 
pattern changes significantly and evolves smoothly from 
that of the ballistic limit [see Fig. ID^a) , which is similar 
to Fig. [IJb)] to that of a classical resistor network [see 
Fig. HJd)]. The evolution is particularly evident when 
plotting the horizontal current in the middle column [as 
shown in Fig. HJc)] for the values of 7 used in Figs. [5a) 
- (d). For 7 = %t^ [corresponding to a dephasing time 
T ~ fi,/(4.9i), see Fig. [DJd)], the current pattern is basi- 
cally indistinguishable from that of the classical resistor 
network (solid purple line). 



In summary, we have shown that nanoscopic quantum 
networks exhibit current eigenmodes with distinct spatial 
current patterns. We demonstrated that the rich variety 
of spatial current patterns arises from the interplay be- 
tween the network's geometry and electronic structure, 
the leads' location and width, and the dephasing time, 
and can be selected via gating or through constrictions. 
These results suggest new venues for custom-designing 
current patterns and their transport properties at the 
nanoscopic, local level. Moreover, we found that a large 
enstrophy of a current pattern generally corresponds to a 
small total current carried by the eigenmode. Finally, we 
demonstrated that with decreasing dephasing time, the 
current patterns evolve smoothly from that of a ballistic 
network to that of a classical resistor network. 
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FIG. 4. Spatial current patterns for (a) 7 = 0.02f^, (b) 7 = 0.18t^, (c) 7 = 0.5t^, and (d) 7 = 8i^ Here, teT = lOOcJo = 
5 X 10~^t, fic = 0, A/i = 0.06i, and leads' Green's function G[ — —iiv. (e) Plot of the horizontal current in the middle column 
of (a) - (d) from left to right and of a classical resistor network (purple line). The current is offset for clarity. 



Appendix A: Dephasing in the presence of phonons 



The Keldysh and retarded Green's function matrices 
are given by 






{f)-'g''{rr+^ 



ph 



G" 



t + s 



ph 



C 



where Ej,^ is the fermionic self-energy matrix arising from 
the electron-phonon interaction. In the limit of temper- 
ature being much larger than the phonon frequency loq 
(i.e., the high-temperature approximation introduced in 
Ref. [l9[), one retains only those terms that contain a 
factor of nsii-^o) and the fermionic self-energy at a site r 
in the network in the self-consistent Born approximation 
is given by 



^r,K 



(w) = ig- 



I^D^MGlfico-.) 



where 



Dq = 2iTT (1 + 2nB(w)) [d{uj + ujq) - 5{lu - ujq)] 

is the Keldysh phonon Green's function, which we assume 
to remain unchanged in the presence of an applied bias, 
and nsii-^) is the Bose distribution function. A further 
simplification is achieved by considering the limit ujq ^- 
in which the self-energy, to leading order in T/uJo, is given 
by 



^K,r 






We next introduce the superoperator D [19] which, 
when operating on a Green's function matrix, returns 
the same matrix with all elements set to zero except for 



the diagonal elements in the network, e.g 



ATi 



[DG 
and thus 



Gj.^, Sr,r' if r lies in the network 
otherwise 



(4) 



S'^'^(cj) = 7DG'-'^ (5) 

We next define the operator U that acts on a matrix X 



UX = CXG" 



(6) 



The solutions of the above Dyson equations are then 
given by 



G^ = U 



l-jDU 



A 



G'' = 



l-gnt + jDG'^ 



g 



(7) 
(8) 



where we defined the diagonal matrix A = gr^9Kga^ ■ 
Note that the only non-zero elements of Arr are those 
where r is a lead site. These elements also contain the 
chemical potentials of the left and right leads. By ex- 
panding the right hand side of the first equation, we ob- 
tain 



G" K \ ^ /Vr 
rr' — 2^ "-^rl 



1 



Aii + 7^QimA„ 



(9) 



m,p 



Ga 
Ir' 



where 



IGfn 



if 1 lies in the network 



otherwise 

Defining next the vector A with A^ == Amm, we finally 



obtain 



rr' — 2^ "-^rl 



1 



I-7Q A 



Ga 



or G = G^YiC^ where the diagonal matrix S is defined 
via 



Sii = 



i-iQ] A 



(10) 
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